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ly 
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Abstract 


This  article  demonstrates  that  the  stress  and  displacement  fields  in 
composites  reinforced  by  hexagonal  fiber  arrays  differ  both  qualitatively  and 
quantitatively  from  results  obtained  for  fiber  and  matrix  layers.  The  layered 
representation  is  a  commonly  used  approximation  for  analyzing  the 
compressive  response  of  composites. 

Finite  element  computational  results  are  presented  for  the  mechanical 
fields  in  uni-axially  reinforced  fibrous  composites  subjected  to  compression 
parallel  to  the  fiber  direction.  The  fibers  are  assumed  to  deflect  in  a  shearing 
mode  deformation  and  both  linear  and  non-linear  responses  are  considered 
for  the  matrix  material. 

The  results  exhibit  a  pronounced  three  dimensional  character,  with 
significant  stress  concentrations  near  the  fiber/matrix  interfaces  and  shear 
stresses  of  magnitudes  equivalent  to  those  of  Xxz-  All  these  features  are 
disregarded  by  the  layered  model  and  may  account  for  its  inadequacy  in 


predicting  compressive  failure. 


Introduction 

With  very  few  exceptions^  the  modelling  of  the  compressive  response 
of  fiber-reinforced  composites  is  based  upon  the  proposition  that  fiber  and 
matrix  regions  can  be  represented  as  layered  domains.  This  simplification 
enabled  Rosen  [1]  to  predict  compressive  failure  by  means  of  a  classical 
buckling  analysis  and  has  been  followed  by  the  great  majority  of  investigators 
since.  Detailed  listings  of  references  can  be  found  in  the  articles  of  Shuart  [2], 
Camponeschi  [3],  Guynn  et  al.  [4]  and  Piggott  [5].  The  above  mentioned 
exceptions,  all  employing  linear  elasticity,  occur  in  the  analyses  of  Sadowsky 
et  al.  [6],  Hermarm  et  al.  [7],  Lanir  and  Fimg  [8]  and  Greszczuk  [9],  where  the 
circular  shape  of  the  fiber  cross-section  was  recognized  in  various  contexts. 
With  the  exception  of  Greszczuk,  the  foregoing  articles  consider  a  single  fiber, 
while  the  complexity  arising  from  the  circular  reinforcement  geometry  is 
simplified  by  Greszczuk  by  means  of  an  assumed  shape  function. 

All  the  above  models  are  deficient  to  various  degrees  in  predicting  the 
compressive  failure  of  fibrous  composites.  Most  predict  failures  at  load  levels 
which  exceed  experimental  observations  and  all  predict  compressive  failure 
loads  which  increase  monotonically  with  fiber  volume  fraction  Vf.  The  latter 
prediction  is  negated  by  experimentally  observed  trends  [10]  and  [11]*. 

The  main  purpose  of  the  present  article  is  to  demonstrate  that  the 
complexities  and  intricacies  of  the  stress  field  in  fiber  reinforced  composites 
subjected  to  compression  cannot  be  represented  by  models  which  employ 
layered  geometries.  Such  models  err  in  two  major  aspects:  (a)  they  overlook 
high  local  stress  concentrations,  especially  near  the  fiber  matrix  interfaces,  and 
(b)  they  account  for  only  one  shear  stress  in  the  matrix  (Xxz)  while  discarding  a 
second  component  (Xyz)  which  is  of  an  equivalent  magnitude.  These 
inadequacies  were  noted  in  a  recent  work  [13],  by  means  of  an  idealized,  linear 
elastic  formulation. 


*  In  a  recent  article  [12],  an  improved  correlation  between  compressive  failure  stress  and  Vf  was 
achieved  by  considering  non-uniform  fiber  spacings.  However,  that  work  still  employed  the 
prevailing  approximation  of  layered  domains. 
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An  impressive  and  conclusive  experimental  evidence  for  the  existence 
of  local  stress  concentrations  and  of  three  dimensional  features  arising  in 
fibrous  composites  under  compression  was  provided  in  a  recent  work  by 
Lankford  [14].  Typical  observations  which  exhibit  the  foregoing  effects  are 
shown  in  Figure  1. 


Basic  Considerations 

Consider  a  composite  material,  uni-directionally  reinforced  by  a  perfect 
hexagonal  fibrous  array  extending  in  the  x-y  plane  as  shown  in  Figure  2a.  Let 
the  fibers  run  parallel  to  the  z  direction,  with  compression  applied  in  the  fiber 
direction.  Due  to  symmetry  it  suffices  to  analyze  the  representative 
rectangular  area  PQRS  shown  in  Figure  2b. 

For  graphite /PEEK  (APC-2)  composites  with  fiber  volume  fraction 
Vf  =  0.4  we  have  a  =  4jj.m  and  b  =  6^,m  in  Figure  2b.  In  addition,  we  presume  a 
lateral  deflection  that  varies  sinusoidally  in  z  with  waviness  length  2L, 
where  L  =  800a.  Such  waviness,  sketched  in  Figure  3a,  is  attributable  to 
processing  non-uniformities,  as  reported  by  several  investigators  [15]  and  [16]. 
Consequently,  it  is  possible  to  confine  the  computations  to  the  rectangular 

block  -  L/2  <  z  <  L/2  with  the  aoss-sectional  area  shown  in  Figure  3b. 

The  computations  assumed  isotropic  fiber  and  employed  values  of 
Ef=  214GPa,  and  Vf  =  0.3.  The  one  dimensional  non-linear  shear  response  of 
the  matrix  shown  by  curve  "a"  in  Figure  4  [15]  was  expressed  by  the  relation 

Y  =  X  /  A  +  (x  /B)^  ,where  A  =  3096,  B  =  169.93,  and  q  =  4.2  (A,  B  and  x  in  MPa). 

The  above  expression  was  extended  to  the  three  dimensional  case  by  means  of 
the  Von-Mises  criterion.  Accordingly,  the  relation  between  the  effective 


strain  £e  and  effective  stress  is  given  by  letting  Be  -  7  and  Cq  =  V^x  in  the 


foregoing  expression.  The  resulting  relation  among  effective  quantities  is 
shown  by  the  curve  "e"  in  Figure  4. 

To  assess  the  significance  of  the  non-linear  matrix  response, 
comparative  computations  were  performed  for  an  assumed  isotropic  linear 
elastic  matrix  behavior  with  Em  =  6.5  GPa  and  Vm  =  0.3. 

In  addition  to  the  foregoing  sets  of  stress-strain  relations  within  the 
fiber  and  the  matrix  regions,  and  the  obvious  requirement  of  stress 
equilibrium  which  are  incorporated  within  the  finite  element  code,  we 
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utilized  the  following  boundary  conditions  for  the  representative  block 
0<x<V3b,  0<y<b, -L/2<z<L/2: 


< 

tl 

o 

N 

II 

o 

II 

o 

on  y  =  0andy  =  b 

u  =  uo  sin  (tcz/L),  v  =  0,  w  =  -  ze^ 

on  X  =  0  and  x  =  b 

(1) 

W  =  i  {L/2)  £2,  Xxz  “  0/  ^yz  ~  0 

on  z  =  ±  L/2 

where  all  symbols  in  equations  (1)  accord  with  standard  notation. 


The  uniform  compressive  strain  £§  in  equations  (1)  is  given  by  ei  = 


A 

L/2' 


where  A  denotes  the  deflection  of  the  plane  z  =  L/2  relative  to  z  =  0,  given  by 


dz  =  ^ua 


(2) 


Upon  selecting  uq  =  8a  the  computations  corresponded  to  a 
compressive  strain 


p  -p0_2li/HQ.V 

41l) 


(3) 


Computational  Details 

In  view  of  the  high  slenderness  ratio,  L/b,  of  the  representative 
rectangular  block  the  computational  scheme  employed  a  coarse  mesh  to  cover 
the  entire  representative  region,  and  a  refined  mesh  to  focus  on  the  region 
-  2a  <  z  <  a. 

The  coarse  mesh  consisted  of  800  layers  normal  to  the  z  direction,  with 
each  layer  further  subdivided  into  16  area  elements  as  shown  in  Figures  5a 
and  5b.  Figure  5a  also  shows  the  deformed  configuration  of  the 
representative  block  with  the  coarse  finite  elements  mesh.  In  the  coarse 
mesh  depicted  in  Figure  5b  each  fiber  is  contained  within  the  three  elements 
at  the  upper  right  and  lower  left  comers. 
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The  stress  and  displacement  fields  were  computed  through  the 
implementation  of  the  commercial  code  ABAQUS  Version  5.2,  using  SUN 
Sparcstations.  The  computations  implemented  the  Von-Mises  plasticity 
model  with  an  assumed  isotropic  hardening.  The  resulting  contours  of  the 
displacement  component  w  at  z  =  0  are  shown  in  Figtire  6,  while  values  of  the 
displacement  components  u  and  v  at  selected  (x,  y)  locations  and  at  z  =  0  are 
listed  in  Table  1.  It  is  obvious  that  the  contours  in  Figure  6  exhibit  strong 
variations  in  both  x  and  y  directions  and  that  the  displacement  v  attains 
amplitudes  comparable  to  those  of  u.  Both  observations  demonstrate 
significant  departures  from  results  predicted  for  layered  geometries. 

To  focus  attention  on  the  stress  and  displacement  fields  at  the  mid¬ 
plane  z  =  0,  where  they  attain  maximal  amplitudes,  the  region  -  2a  <  z  <  a  was 
subdivided  by  a  finer  mesh  into  30  layers  normal  to  z,  with  each  layer 
consisting  of  116  elements  as  shown  in  Figure  7.  The  refined  mesh  allots  each 
fiber  39  elements  adjacent  to  the  upper  right  and  lower  left  corners  in  that 
figure. 

The  asymmetry  in  the  range  of  z  about  z  =  0  in  the  refined  mesh  was 
chosen  deliberately  to  assess  the  adequacy  of  the  selected  range  -  2a  <  z  <  a. 
This  assessment  is  provided  by  examining  the  magnitudes  of  relative 
deviations  in  computed  values  from  the  obvious  requirements  of  symmetry 
about  rays  emanating  from  the  central  point  x  =  V3  b/2,  y  =  b/2.  (The  above 
range  turned  out  to  be  adequate  since  the  computational  results  met  those 
requirements  to  a  high  degree  of  accuracy). 

Accordingly,  the  displacements  computed  with  the  coarse  mesh  at 
z  =  -  2a  and  z  =  a,  as  well  as  their  interpolated  values  at  locations  x,  y  which 
matched  positions  on  the  refined  mesh,  were  employed  as  prescribed 
boundary  values  for  the  detailed  computation.  For  purpose  of  illustration, 
the  boimdary  values  of  w  (x,  y)  at  z  =  a  are  sketched  in  Figure  8. 

Comparisons  with  predictions  derived  for  layered  geometries  were 
obtained  by  considering  a  rectangular  block  -  L/2<z<L/2  with  an  inner 
"matrix  layer"  and  outer  "fiber  layers"  contained  within  0  <  x  <  t/2, 
fSh-  \i/2  <x<f3h,  0<y<b,  respectively,  as  shown  in  Figure  9.  The  equality 

2 

of  Vf  requires  that  t  =  .  The  deformed  shape  of  the  layered  geometry  is 

2b 

shown  in  Figure  10. 
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Results  and  Discussion 


Contours  of  the  displacement  w  and  the  stresses  '^yz  ®z 
shown  in  Figures  11  and  12  for  the  linear  case  and  in  Figure  13  for  the  non¬ 
linear  matrix  response  depicted  in  Figure  4. 

Figures  11  and  13  exhibit  contoms  of  the  above  mentioned  components 
at  the  mid-plane  z  =  0  while  Figure  12  provides  perspective  views  of  the 
deformed  blocks  between  z  =  -  a  and  z  =  a/2.  For  additional  clarity  the 
displacements  w(x,  0 , 0)  and  w(x,  b,  0)  are  plotted  vs.  (o  <  x  <  b)  in  Figures 
14a  and  14b. 

The  results  shown  in  Figures  11  through  14  are  contrasted  against 
predictions  which  correspond  to  layered  "fiber"  and  "matrix"  regions  shown 
in  Figures  15  and  16.  Comparative  values  of  the  largest  stress  components  are 
given  in  Table  2.  That  table  lists  values  for  circular  and  layered  geometries, 
denoted  by  rows  labeled  "C"  and  "L",  respectively,  and  contains  comparisons 
between  linear  and  plastic  matrix  response. 

The  predominant  conclusions  to  be  drawn  from  the  current  results  are: 

(1)  It  is  erroneous  to  represent  fibrous  arrays  of  fiber  reinforced  materials 
by  equivalent  layered  media  for  evaluating  compressive  response. 

(2)  The  stress  and  displacement  fields  in  fibrous  arrays  subjected  to 
compression  are  markedly  non-uniform  and  exhibit  highly  localized  effects, 
especially  near  the  fiber/ matrix  interfaces. 

(3)  The  shear  stress  Xy^,  which  is  disregarded  in  the  layered  representation, 
attains  magnitudes  which  are  equivalent  to  those  of 

(4)  Compressive  failure  of  fiber  reinforced  composites  may  arise  due  to 
localized  effects  which  are  overlooked  by  the  layered  model. 
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Fig.  2a.  Hexagonally  arrayed  fibers. 


Fig.  2b.  The  representative  cross-section  of  the  hexagonal  array. 


Fig.  3a.  A  sketch  of  fiber  waviness. 


Fig.  3b.  The  representative  volumetric  block. 
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Fig.  5a.  Deformation  overview  and  coarse  mesh  geometry. 
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Fig.  8.  A  sketch  of  the  distributed  nodal  displacements,  computed 
with  the  coarse  model,  applied  as  boundary- 
conditions  of  the  top  surface  of  the  fine  mesh. 
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Fig.  9.  The  equivalent  layered  fiber/matrix  geometry. 


Fig.  10.  The  deformed  shape  of  the  layered  geometry  under  analogous 

compressive  circumstances. 


17 


(a)  The  displacement  w 
Contours  "1"/  "2",  and  "3"  correspond 
to  w  =  -  9.65  X  10‘8,  8.35  x  lO'^,  and 
9.57  X  10'^  m,  respectively. 


=  -  27.7, 1.35,  and  10.3  MPa,  respectively. 


=  35.4,  209,  and  248  MPa,  respectively. 


(d)  Shear  stress  Xy^ 

Contours  "1",  ”2",  and  "3"  correspond  to 
Xyz  =  -  2.11,  7.3,  and  101  MPa,  respectively. 


Fig.  11.  Linear  results:  Displacement  and  stress  contours  in  the  mid 
plane  (z  =  0)  with  circular  fibers  geometry. 
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(a)  The  displacement  w 
Contours  "1",  "2",  and  "3"  correspond 
to  w  =  -  1  X  10'7  9.43  X  10'9,  and 
lx  10'^  m,  respectively. 


to  (522  -  ■  ~  ^5.9, 93.3,  and 

147  MPa,  respectively. 


to  =  24.6,  55.0, 100.0,  and 
19.2  MPa,  respectively. 


(d)  Shear  stress  Xyz 
Contours  "1",  "2",  "3",  and  "4”  correspond 
to  Xyz  =  -  35.8,  -  6.1, 5.82,  and 
29.6  MPa,  respectively. 


Fig.  13.  Nonlinear  results:  Displacement  and  stress  contours  in  the  mid 
plane  (z  =  0)  with  circular  fiber  geometry. 
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(a)  The  displacement  w 
Contours  "1",  "2",  and  "3"  correspond 
to  w  =  -  5.3  X  lO'S,  4.7  X  10^,  and 
5.28  X  lO'S  m,  respectively. 


(b)  Normal  stress 
Contour  "1"  corresponds  to 
=  -  69.9  MPa. 


(c)  Shear  stress 

Contours  "1",  "2",  and  "3"  correspond  to 
=  128, 129,  and  130  MPa. 


Fig.  15.  Linear  results:  Displacement  and  stress  contours  in  the 
mid  plane  of  the  layered  geometry. 

[Note:  The  dashed  lines  represent  the  interfaces  between  "fiber”  and  "matrix”  layers.] 
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(a)  The  displacement  w 
Contours  "1",  "2",  and  "3"  correspond 
to  w  =  -  5.44  X  10*3,  5.0  X  and 
5.42  X  10*3  m,  respectively. 


(b)  Normal  stress  a^z 
Contours  "1",  "2",  and  "3"  correspond  to 
G^z  =  -  103,  -  60.2,  and  8.82  MPa,  respectively. 


(c)  Shear  stress 

Contours  ”1"  and  ”2"  correspond  to 
=  52.7  and  77.8  MPa,  respectively. 


Fie.  16.  Nonlinear  results;  Displacement  and  stress  contours  m 
the  mid  plane  of  the  layered  geometry. 


[Note;  The 


dashed  lines  represent  the  interfaces  between  "fiber"  and  "matrix"  layers.] 
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Table  1:  Computational  results  for  the  displacement  components  u,  v  and  w 
(in  meters)  at  various  positions  (x;  y)  in  the  mid  plane  z  =  0.  • 
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Table  2:  A  comparative  listing  of  maximal  and  minimal  stresses  at  the 
mid-plane  z  =  0,  for  circular  fibrous  (C)  and  layered  (L)  geometries. 
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